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Abstract 



Two general methods for establishing the logarithmic behavior of recursively defined sequences of 
real numbers are presented. One is the interlacing method, and the other one is based on calculus. 
Both methods are used to prove logarithmic behavior of some combinatorially relevant sequences, such 
as Motzkin and Schroder numbers, sequences of values of some classic orthogonal polynomials, and 
many others. The calculus method extends also to two- (or more- ) indexed sequences. 

Keywords: log-concavity, log-convexity, special combinatorial numbers, calculus, integer sequences 

AMS subject classifications: 05A20, 05A16, 05E35, 11B83, 11B37, 11B39, 11B68 



2 



1 Introduction 



Let a(n), n > 0, be a sequence of positive real numbers. We want to examine the rate of growth of 
this sequence, i.e. to examine whether the quotient ^^^"''^^ decreases, increases or remains constant. 
In other words, we want to see whether the sequence is log-concave, i.e. a(n)^ > a{n — l)a(n + 1), 
log-convex, i.e. a(n)^ < a{n — l)a{n + 1), or log-straight (or geometric), i.e. a(n)^ = a{n — 
l)a{n + 1) for all n > 1. Under log-behavior we also sometimes include log-Fibonacci behavior, 
meaning sign[a{n)^ — a{n — l)a(n + 1)] = sign{—l)^ (or (—1)"+-'^). It is of great interest, especially 
in combinatorics, as it can be seen from many examples in j25j . to know the log-behavior of a given 
sequence. It is, in fact, just one instance of the whole paradigm of "positivity questions" ( 28 ). 

If a(n) has a combinatorial meaning it would be ideally to provide a combinatorial proof of its log- 
behavior. For example, if we want to prove that a{n) is log-convex and if we know that a{n) = \S{n)\, 
where S{n) is a certain finite set, then we would like to find an injection S{n) x S{n) S{n — 1) x 
S{n + 1), or a surjection S{n — 1) x S{n -|- 1) ^ S{n) x S{n), and similarly for log-concavity. It is 
usually a hard task to find such a (natural) injection or surjection. Still, examples of this type include 
binomial coefficients, Motzkin numbers ([HI) and permutations with a prescribed number of runs ([Sj). 
Of course, the explicit formulae give another possibility to prove results of this type, but they are 
rarely on disposal. Instead, other methods for proving such inequalities have been developed, e.g. see 

m, m, m or u- 

In this paper besides using old methods to prove some new results on log-behavior, we shall also 
introduce some new methods and use them to prove log-behavior of certain interesting combinatorial 
sequences, and apply this method to other sequences, the most prominent example being values of 
classical orthogonal polynomials. 
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2 Log-behavior of some sequences using known results 



Let us quote some known results and apply them to examine the log-behavior of certain combinatorial 
and other sequences. 

Lemma 2.1 (Newton's lemma) 

Let P{x) = Ylk=o '^kX^ be a real polynomial whose all roots are real numbers. Then the coefficients 
of P{x) are log-concave, i.e. a| > afc_iafe+i, A; = 1, . . . , n — 1. Moreover, the (finite) sequence is 

\k) 

log-concave in A;. ■ 

Let us briefly recall how to apply this lemma to binomial coefficients and Stirling numbers c(n, k) of 
the first kind (the number of permutations on the set [n] = {1,2, ... ,n} with exactly k cycles) and 
Stirling numbers S{n,k) of the second kind (the number of partitions of [n] into exactly k blocks). 
The following formulae are well known: 



k=0 

n 

= ^ c(n, A;)x^ (2.2) 
fe=o 

n 

x" = ^ S{n, k)x^, (2.3) 



where x- = x{x — 1) . . . (a; — fc + 1) is the fc-th falling power, and x^ = x{x + 1) . . . {x + k — 1) the 
k-th rising power of x. Prom (2.1) and (2.2) we see that {x + 1)"^ and x"^ have only real roots. So, 
by Newton's lemma we conclude that the sequences (^) and c{n,k) are log-concave. The case of the 
sequence S{n, k) is a bit more involved. Let 

n 

Pn{x) = Y,S{n,k)xK 

k=0 
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From Po{x) = 1 and from the basic recursion 

S{n, k) =S{n-l,k-l) + kS{n - 1, A:), 

it follows at once that 

Pn{x) = X[P^{X) + Pn-l{x)]. 

The function Qn{x) = Pn{x)e^ has the same roots as Pnix) and it is easy to verify Qn{x) = xQ'j^{x). 
By induction on n and using Rolle's theorem it follows easily that Qn, and hence P„, have only real 
and non-positive roots. So, we conclude: 

Theorem 2.2 

The sequences {X)k>Q^ ('^(^'^))fc>o i^i^^ ^))k>o log-concave. Hence, they are unimodal. ■ 
An inductive proof of Theorem 2.2 is given in |22j . 

The next easy lemma is sometimes useful in proving log-convexity results. 
Lemma 2.3 

Let / : [a, 6] — > M be a positive, continuous function, and 

I„ = f f{xTdx, n>l. 

J a 

Then (/n)„>2 is a log-convex sequence. 
Proof 

By Cauchy-Schwarz inequality, we have 

■ 
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As an example, we apply this lemma to Legendre polynomials Pn{x). It is well known (e.g. j31j ) that 
the following Laplace formula holds: 

1 

P„(x) = - / (x + Va;2 - lcos(p)''dip. (2.4) 
71" Jo 

Hence, from Lemma 2.3 and (2.4) we obtain 
Theorem 2.4 

The values Pn{x), n > 0, for x > 1 are log-convex. ■ 
Another proof of this fact will be presented in Section 4. 

We say that a sequence (an)„>o has no internal zeros if there do not exist integers < i < j < k such 
that Oj 7^ 0, aj = 0, Ofc 7^ 0. 

Theorem 2.5 (Bender-Canfield, see [5]) 

Let 1, ai, a2, . . . be a log-concave sequence of nonnegative real numbers with no internal zeros and let 
(^n)n>o sequence defined by 

= exp^afc— . (2.5) 

n>0 \A.>1 / 

Then the sequence (6n)n>o is log-convex and (7^)„>o log-concave. ■ 

As our first application of this theorem consider the Bell numbers {Bn)n>o- Bq = 1, and Bn is the 
number of partitions of an n-set. It is well known that the exponential generating function for (-Bn)„>o 
is given by 

/ 

^ ^"7^ " ^^^^^"^ - 1) = ea;p 5Z 

n>0 ■ \A;>1 

By taking = -^^^Yyii k>\, and checking that the sequence 1, oi, 02, . . . is log-concave, we conclude 
from Theorem 2.5 that {Bn)n>o log-convex and (■^)„>g is log-concave sequence. More generally, 
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for an integer I > 2, define expi to be I times iterated exponential function, i.e. 

expi{x) = exp{exp{. . . {exp{x)) . . .), 
exp written I times. Define the sequence (bn^) by 

V /n>0 



n>0 

and Bell numbers of order I by 

r(0 - (2 6) 

" ~expi{Oy ^'-^^ 

(2) 

So, Bn = Bn are ordinary Bell numbers. Now it is not hard to prove by induction on I the following 
result. 

Theorem 2.6 

For any fixed I > 2, the Bell numbers of order Z, i.e. the sequence [Bn^j is log-convex, and the 

V / n>0 

sequence is log-concave. ■ 

The following lemma is an easy consequence of the definition of log-convexity and log-concavity. 
Lemma 2.7 

Suppose (an)„>o is a positive, log-convex sequence and ao = 1. Then an«m < «n+m- If) in addition, 
(^)n>o log-concave, then Qn+m < ■ 

Now Theorem 2.6 and Lemma 2.7 imply semi-additivity inequalities for Bell numbers of order I. 

Corollary 2.8 



Bll'BiJ}<BZn<{^^^")Bll^B(^, m,n>0. 
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The next apphcation of the Bender- Canfield theorem concerns the number of certain permutations, 
i.e. elements vr G the symmetric group of permutations on n letters. For a fixed integer k > 1 
let Cfc(n) be the number of all permutations from Sn that have cycles of length not exceeding k. By 
definition, Cfc(O) = 1. 

Theorem 2.9 

The sequence (cfc(n))^>Q is log-convex and the sequence (j^^:^^ is log-concave. 
Proof 

It is well known (see, e.g. that the exponential generating function Ck{x) of the sequence 

(cfc(™))„>0 is given by 

Hence the sequence {ak)k>o fro™ the Bender-Canfield theorem in (2.5) is as follows: oq = l,ai = 
l,...,afc = l,ak+i = 0, aA:+2 = 0,..., which has no internal zeros and is obviously a log-concave 
sequence. Hence the claim follows from Theorem 2.6. ■ 

The requirement that the sequence (a^) from Bender-Canfield theorem does not have internal zeros is 
essential, and in general can not be weakened. As an illustration, let us consider a class of sequences 
related to Cfc(n). Again, for a fixed integer > 1 let efc(n) be the number of all permutations vr from 
S„ such that tt'^ = id, for n > 1 and 6^(0) = 1. Then the exponential generating function Ek{x) of 
this sequence is given by (see, e.g. j^T]) 

Eu{x)=exp{Y,-y (2.7) 

and the corresponding recurrence is 



efc(n) = ^(n - iy—Lek{n - j) 
j\k 



with appropriate initial conditions. The sequence (a„) from (2.5) is given by oq = ai = 1, aj = 1 
if j divides k and aj = otherwise. This binary sequence 1, 1, 02, . . . , a^, 0, 0, . . . is log-concave if 
and only if it does not contain 1,0, 1 as a subsequence, but it contains internal zeros for all k > 2. 
Taking, for example, = 5, the first few terms of 65(71) being 1,1,1,1,1,25,145,505,1345,..., we 
can easily see that this sequence does not exhibit any logarithmically definite behavior, although the 
sequence 1,1,0,0,0,1,0,... is log-concave. For higher values of k, sequences eyt(n) log-behave even 
more chaotically. 



3 The interlacing (or "sandwich") method. Secondary structures. 

Let a{n), n > 0, be a sequence of positive numbers defined by a homogeneous linear recurrence, say 
of second order: 

a(n) = R{n)a{n - 1) + S{n)a{n - 2), (3.1) 
where R and S are known functions, together with given initial values a(0) = Oq, and a(l) = ai. 

Let our task be to examine the rate of growth of a(n). We define the sequence of consecutive quotients 

q{n)= / n>l. 
a{n — 1) 

Dividing (3.1) by a{n — 1) we obtain the recurrence 

^H = ^H + ^^31)' (3-2) 

with initial condition q{l) = ^ '■= bi. The log-concavity or log-convexity of (a(n)) is equivalent, 
respectively, to q{n) > q{n + 1) or q{n) < q{n + 1), for all n > 1. So, what we want to see is whether 
the sequence {q{n))^~^i decreases or increases. To prove that {q{n)) increases, it is enough to find an 
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increasing sequence (6(n)) such that 

b{n) < q{n) < b{n + 1) (3.3) 

holds for all n > 1, or at least for all n > no for some uq. Then we can conclude that (a(n)) behaves 
log-convex at least from some place on. Analogously for log-concavity. This "sandwich method" or 
"interlacing method" works in some simple cases, but often it is very hard to hit the right sequence 
{b{n)) which is simple enough. In the rest of this section we show how this method works for some 
combinatorially important sequences. We also show some consequences of the obtained results. 

Example 3.1 (Derangements) 

Let Dn be the number of derangements on n objects, i.e. the number of permutations tt G S„ without 
fixed points, for > 1, and Dq := 1. It is well known (and easy to prove) that the following recurrence 
holds: 

Dn = {n - l)[Dn-i + Dn-2], (3.4) 

with initial conditions Dq = 1, Di = 0. Then D2 = I, D3 = 2, = 9, D5 = 44, Dq = 265 etc., and 
we expect from these initial values that < for n> S. Indeed, divide (3.4) by -Dn-i S'ld 

denote q{n) = -^^^ Then 



q{n) = (n - 1) 



1 

1 + 



(3.5) 



q{n — 1) ' 

with q{3) = 2, ^^(4) = 9/2. Let b{n) = n — 1/2. It is easy to check by induction on n and using (3.5) 
that 

6(n) < q{n) < 6(n + 1), 

for all n > 4. Since {b{n)) is clearly increasing and q{3) < g(4), we conclude that (D„) is log-convex 
for n > 3. ■ 
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Example 3.2 

Let T2 (n) denotes the number ofnxn symmetric No-matrices with every row (and hence every column) 
sum equal to 2 with trace zero (i.e. all main-diagonal entries are zero) (Example 5.2.8 in 27^). The 
exponential generating function of T2(n) is given by 

1 / x'^ X \ 

^^") = 7r=^^"HT-2j- 

The numbers 72 (n) satisfy the recurrence 

T2(n) = (n - l)r2(n - 1) + (n - l)T2(n -2)-(^''~^^ T2{n - 3) 

with the initial conditions T2(0) = 1, T2(l) = 0, T2(2) = T2(3) = 1. The corresponding recurrence for 
successive quotients q2in) = T2{n)/T2{n — 1) is given by 

n — 1 f n — 1\ 1 

q2{n) = (n - 1) + — — - — -— — , n > 5, 

q2{n-l) V 2 J q2{n - l)q2{n - 2) 

with the initial conditions (72(3) = 1, 92(4) = 6. Tabulating the first few values of q2{n), we see 
that, after some initial fluctuations, this sequence seems to behave like n — ^. So, we guess that 
n — 1 < q2in) < n, and indeed, this follows easily by induction on n for n > 6. Hence, the sequence 
(72(n))„>g is log-convex. ■ 

For our next application we need some preparations. A Motzkin path of length n is a lattice path 
in (x,y)-plane from (0,0) to (n, 0) with steps (1,1) (or Up), (1,-1) (or Down) and (1,0) (or Level), 
never falling below the x-axis. Denote by M{n) the set of all Motzkin paths of length n. The number 
Mn = |A^(n)| is the n-th Motzkin number. By definition, Mq = 1. 

A handful of other combinatorial interpretations of M„ are listed in Ex. 6.38 in 27 . A typical member 
of the Motzkin family A^(20) is shown in Fig. 1 as a "landscape path". A peak of a Motzkin path 
is a place where an Up step is immediately followed by a Down step. A plateau of length I is a 

11 




peak 
terrace \ 

/ / \ trench 



plateau 



valley 



plain 



\' A 



(0,0) 



(20,0) 



Figure 1: A Motzkin path of length 20 

sequence of I consecutive Level steps, immediately preceded by an Up step and immediately followed 
by a Down step. Similarly we define a terrace, trench, valley and plain of a Motzkin path. 

A Motzkin path without any Level steps is called a Dyek path or a "mountain path". It is well 
known that the set P(n) of all Dyck paths of length 2n is enumerated by Catalan numbers C„, 
i.e. |X'(n)| = Cn = i^n)- I^o™ the explicit formula it follows immediately that < Cn~iCn+i, 
i.e. Catalan numbers are log-convex. It is also easy to find a simple combinatorial proof of this fact. 
By counting the number of Level steps on a Motzkin path, it can be easily shown that Catalan and 
Motzkin numbers are related as follows: 



Proposition 3.2 

a) The Motzkin numbers satisfy the following convolutive recursion: 

n-l 
k=0 

b) The generating function of (M„) is given by 





(3.7) 
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c) M„'s satisfy the short recursion 

(n + 2)Mn = (2n + l)Mn-i + 3(n - l)Mn-2- (3.8) 

Proof 

a) Mq = 1 and Mi = 1. A Motzkin path of length n + 1 either starts by a Level step and then can 
proceed in M„ ways, or starts by an Up step and returns for the first time to the x-axis after k + 1 
steps (to the point {k + 2,0)). The number of latter is equal to the number of pairs {Pi,P2), where 
Pi is a (translated) Motzkin path from (1, 1) to {k + 1, 1) which is not below the line y = 1, and P2 
is a (translated) Motzkin path from (k + 2, 0) to (n + 1, 0). The number of paths Pi is equal to M^, 
while the number of paths P2 is equal to the number of Motzkin paths onn + 1 — (A: + 2) = n — k — 1 
steps, and this is Mn-k-i- Thus a) follows. 

b) If we multiply (3.6) by x""^^ and sum over n > 0, we get the functional equation a;^M^(x) + (x — 
l)M(x) + 1 = 0, and since M(0) = Mq = 1, we obtain (3.7). 

c) The generating function M{x) is algebraic, hence D-finite. Therefore (M„) is P-recursive. Now 
from Eq(6.38) in '27 for the polynomial A(x) = 1 — 2x — 3x^ and r = 2, ao = 1, ai = —2, 02 = —3, 
c? = 2, we get the claim. ■ 

Theorem 3.3 

a) The sequence M„ of Motzkin numbers is log-convex. 

b) Mn/Mn-i < 3, for all n > 1. 

c) The sequence M„/M„_i is convergent and lim,„^oo -^n/-^n-i = 3. 
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Proof 

a) Divide (3.8) by (n + 2)M„_i and let q{n) := M„/M„_i. Then we obtain 



g(n) = 



n + 2 



3(n 
2n + l+ ^ 



g(n- 1)_ 



n > 2, (3.9) 



with the initial condition q{l) = 1. Then q{2) = 2. Now define the sequence 6„ = S ^^^g . This 
sequence is obviously increasing. We claim that bn < q{n) < bn+i, for all n > 3. This is obviously 
true for n = 3. Let n > 3. By inductive hypothesis, we have 

, , 2n + 1 3(n - 1) 1 ^ 2n + 1 3(n - 1) 1 

(7(77,) = ~\~ — ^ ^ ~\~ — ^ 

n + 2 n + 2 q{n-l) ~ n + 2 n + 2 bn 

But 

2n + 1 3(n - 1) 1 ^ _ ^ " - 3 ^ ^ 

n + 2"^ n + 2 " ~ 2 n(n + 2)(2n + 3) " ' 

for all n > 3, and the inequality q{n) > bn follows. On the other hand, by the inductive hypothesis 
bn-i < (l{n — 1), we have 

, , 2n+l 3(n-l) 1 
q{n) < — — + 



n + 2 n + 2 bn-i 
Subtracting bn+i from the right hand side, we get 

2n + 1 3(n - 1) 1 _ _ _ 1 1 
n + 2 ^ n + 2 bn-i ~ ~ ~ 2 (n + 2)(2n + 5) " ' 

and the inequality q{n) < follows. So, the claim is proved by induction. Hence, the sequence 
Mn/Mn-i is increasing, and the Motzkin sequence is log-convex. 

b) We have proved in a) that M„/M„_i < 5„ = < 3. 

c) By a) and b) it follows that q{n) = M„/M„_i is an increasing and bounded sequence: 2 < q{n) < 3, 
hence convergent. Passing to limit in (3.9), or passing to limit in the sandwich inequality bn < q{n) < 
bn+i above, the last claim in c) follows. ■ 
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Corollary 3.4 

a) The sequence {Mn/n\) is log-concave. 

b) M2 < Mn-iMn+i < (1 + ^) Ml for all n>l. 

c) M^Mn < M^+n < C"^") M^Mn, for all m,n > 0. 

Proof 

a) The log-concavity of Mn/n\ is equivalent to ^j^"*"^ < jvf^"^ • So, we have to prove that q{n + l) < 
^^q{n). We know that the sequence q{n) is increasing. Starting from the short recursion (3.9) for 
q{n + 1), we have 

2n + 3 3n 1 2n + 3 3n 1 



n + 3 n + Sq{n) n + S n + 3q{n — l) 
n + 2 n 



n + 3 n — 1 



{2n + 3){n-l)2n + l _^3{n-l) 1 



2n{n + l) n + 2 n + 2 q{n-l)_' 

The claim now follows by noting that the term ^^^^(n+i) clearly less than one for all n > 1, and 
that the inequality ^^^^ < is valid for all n > 3. The validity of our claim can easily be checked 
for 1 < n < 3. 

b) This follows from log-convexity of Motzkin numbers and a). 

c) A simple combinatorial proof of the left inequality follows from the fact that a concatenation of 
two Motzkin paths of lengths m and n, respectively, is again a valid Motzkin path of length m + n. 
To prove the right inequality, start from q{n) > ^^(?(n + 1). Using this inequality repeatedly, we get 



Mo-2Mi-3M2~"'~m + n Mm+n-i ' 



for all ra > 0, m > 1. 
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Hence, for any < j < m — 1, we have 

Mj+i ^ i + 1 Mm+n 

Mj ~ m + n Mm+n~i' 

From this we get 

M1M2 Mm ^ f 1 Mn+l \ f 2 Mn+2 \ f m Mm+n \ 

Mo Ml'" Mm-i -\n+l Mn )\n + 2 M„+J " \m + n M„+„_J ' 
and after the cancellation: 

Mm ^ m!n! Mm+n 
Mq - (m + n)! M„ 

Since Mq = 1, we get the claim. The case m = 0, n > is trivial. ■ 
Remark 3.5 

An algebraic proof of log-convexity of Motzkin numbers was given in 1 , and a combinatorial proof 
in [S]. We shall give yet another ("calculus") proof later in Section 4. 

Now we proceed in applying the "sandwich method" to combinatorial structures that generalize both 
Dyck and Motzkin structures in a sense that they are counted by Catalan and Motzkin numbers in 
some special cases. They are called secondary structures and come from molecular biology (see, e.g. 
UHl, IS], El, CZl, m, CI)- More details are in fE] and HI]. 

Let n and I are integers, n > 1, / > 0. A secondary structure of size n and rank / is a labeled 
graph S on the vertex set V{S) = [n] = {1,2, ... ,n} whose edge set E{S) consists of two disjoint 
subsets, P{S) and H{S), satisfying the following conditions: 

(a) {i,i + 1} G P{S), for ah 1 < i < n - 1; 

(b) {i,j} G H{S) and {i,k} e H{S) ^ j = k; 

(c) {i,j} G H{S)^\i-j\>l; 

(d) {i,j} £ H{S),{p,q} £ H{S) andi<p<j^i<q<j. 
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The vertices of a secondary structure are (in biology) usually called bases, the edges in P{S) are called 
p-bonds and those in H{S) /i-bonds. A secondary structure S with H(S) = is called trivial. The 
number of /t-bonds of S is called its order. An example of a secondary structure of size 12, rank 2 
and order 3 is shown below in Fig. 2. Note that every /i-bond "leaps" over at least two bases. We 
p-bonds/'"'^ h-bonds 




1 2 3 4 5 6 7 8 9 10 11 12 

Figure 2: An example of a secondary structure 

denote the set of all secondary structures of size n and rank I by S^^\n), and the set of such structures 
of order k hyS^\n). The respective cardinalities we denote by S^''^{n) and S^\n). By definition, we 
put S(')(0) = 1, for all I. 

Another interpretation of secondary structures is as follows. Denote by M.^^\n) the set of all Motzkin 
paths in M.{n) whose every plateau is at least I steps long. For I = 0, M.^'^\n) = M.(n), i.e. all 
Motzkin paths on n steps. 

Proposition 3.6 

There is a bijection between S^''\n) and M.^''\n) for all n > 1, Z > 0. 



Proof 

We shall only briefly describe the correspondence 5(')(n) —>■ M.^''\n). Take a secondary structure 
S G S^^\n) and scan it from left to right. To each base, starting from 1, we assign a step in a lattice 
path as follows. To any unpaired base (i.e. base in which no /i-bond starts or ends) we assign a Level 
step. To a base in which a new /i-bond starts assign an Up step, and to each base in which an already 
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encountered /i-bond terminates, assign a Down step. The obtained path P is in M.^^\n) and it is not 
hard to check that S i — > P is a bijection. ■ 

An example of the correspondence (S(^)(10) < — > A^(^)(10) is shown in Fig. 3. Therefore, the rank 




Figure 3: A secondary structure and the corresponding Motzkin path 

secondary structures, i.e the border case, are just all Motzkin paths, hence S^^\n) = Note that 
for Z > 1 a secondary structure of rank I is a simple graph. 

If we allow rank to degenerate to the value of Z = — 1, then an /i-bond can terminate in the very vertex 
it starts from, i.e. we allow loops. 

Proposition 3.7 

There is a bijection between the set S^~^\n) of all secondary structures of size n and rank —1 and 
the set P(n + 1) of all Dyck paths on 2(n + 1) steps. Hence, S^~^\n) = Cn+i, for all n > 0. 

Proof 

Again, we shall only briefly describe how to assign a member of S^"^^ (n) to a Dyck path P G V{n + 1). 
So, take a Dyck path on 2(n + 1) steps. Discard the first and the last step and divide the remaining 
path in groups of two consecutive steps. Assign to each group a base in a secondary structure according 
to the following rule. To a group of two Up steps assign a base in which an /i-bond starts, to a group 
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of two Down steps assign a base in which the last started /i-bond terminates. To a group {Up, Down) 
assign a base with a loop attached to it, and, finally, to a group (Down, Up) assign an unpaired base. 



Figure 4: The correspondence between secondary structures of rank -1 and Dyck paths 

interpretation of secondary structures is related to pattern avoidance in permutations. A pattern is 
a permutation cr G S^, and a permutation vr G S„ avoids a if there is no subsequence in tt whose 
members are in the same relative order as the members of a. It is well known ( JHl) that the number 
of permutations from E„ avoiding o" G S3 is equal to the Catalan number C„, for all patterns a G E3. 

The concept of pattern avoidance was generalized in [S], by allowing the requirement that two letters 
adjacent in a pattern must be adjacent in the permutation. An example of a generalized pattern is 
1 — 32, where an 1 — 32 sub- word of a permutation vr = 0102 ... a„ is any sub- word ajOjaj+i such that 
i < j and < aj+i < aj. Generalized pattern avoidance is treated in more detail in jnj, where it 
is shown that the permutations from S„ that avoid both 1 — 23 and 13 — 2 are enumerated by the 
Motzkin numbers. 

Proposition 3.8 

For any / > 0, there is a bijection between S^^\n) and the set of all permutations from S„ that avoid 




An example of the correspondence I'(8) < — > ^\'^) is shown in Fig. 4. The third combinatorial 



1 
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{1 - 23, 13 - 2, ij}, where j>i + l. 



Proof 

This bijection is just the bijection given in the proof of Proposition 24 in 'W\- It is easy to see that 
any pattern i,i + k in a permutation a £ $]„ avoiding {1 — 23, 13 — 2} generates a plateau of length 
A; — 1 in the corresponding Motzkin path. So, the claim follows from Proposition 3.6. ■ 

Let us find now the recurrences and generating function for the general secondary structure numbers 
analogous to those for the Motzkin numbers (the border case) as was shown in Proposition 3.2. 

Proposition 3.9 

a) For any fixed integer I > —1, the numbers S^''\n) of secondary structures of rank I satisfy the 
following convolutive recurrence 

n-l 

5(')(n+ 1) = 5(')(n) + ^5(')(m)5(')(n-m- 1), n>l + l, (3.10) 

m=l 

together with the initial conditions 

5(0(0) = 5(0(1) = . . . = 5(0(; + 1) = 1. (3.11) 

b) The generating function 5(')(x) = X^„>o ^('^(n)^" is given by 



-uJi{x) — ^ u>f{x) — 4x2 



where 



Si{x) = , (3.12) 



uji[x) = x-l-x'^ - x'+^ = 2x - (1 + X + x^ + . . . + x'+^) = 1 - x^- — —. (3.13) 

1 — X 



c) The sequence ('S'^'H"'))„>o satisfies the following "short" recursion 
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+ 2)5«(n) = ^ A«(n,fc)5«(n- A:), (3.14) 



[n 

k=l 
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where 

-i(A;-3)(2n + 4-3/c) , l<k<l + l 

A^^\n,k) = <. -i(/-3)(2n-3/-2) , k = l + 2 

^ -i(2/ + 3-/c)(2n + 4-3A:) , / + 3 < /c < 2/ + 2. 

Proof 

a) Clearly, for n < / + 1 there is only one (the trivial) secondary structure of rank I and size n, and 
hence initial conditions (3.11) hold. Let n > / + 1. A secondary structure on n + 1 bases either does 
not contain an /i-bond starting at the base 1, in which case there are S^^\n) such structures, or it 
has an /i-bond from the base 1 to some base m + 2, at least / bases apart. In this case, there are 
5'(0(j72)5'(0(7j — — I'j such structures, and (3.10) follows. 

b) Multiplying (3.10) by x""''^, summing over n > and taking into account the initial conditions 
(3.11), we obtain the functional equation 

x^[Siix)f +ujiix)Si{x) + l = 0, 

and this, in turn, implies (3.12). 

c) As in the proof of Proposition 3.2 c), we use D-finiteness of the above generating function and again 
formula (6.38) in ;27] with d = 2, r = 2/ + 2 to obtain the claim. More details are in [llj and ^2]- ^ 

Let us write down explicitly the recurrences (3.14) in cases / = 1, 2, 3 (of course, for / = 0, it coincides 
with (3.8)): 

5«(n) = ^5«(n - 1) + ^5«(n - 2) + ^5«(n - 3) - ^^^(n - 4), (3.15) 
n+2 n+1 n+2 n+1 

5^(0) = S«(l) = 5^(2) = 1,5(1) (3) = 2; 
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(3.16) 

5(2) (0) = (1) = 5(2) (2) = 5(2) (3) = 1, 5(2) (4) = 2, 5(2) (5) = 4. 



5(3) (^) = ^!i±i5(3)(n - 1) + !^5(3)(n - 2) 



n + 2 
3n-21 



n + 2 



5(3)(n-6) 



n + 2 
2n- 17 



n - 4 
n + 2 



5(3) (n- 4) 



n + 2 



5(^)(n-7)-^5(3)(n-8), 



(3.17) 



5(3) (0) = 5(3) (1) = 5(3) (2) = 5(3) (3) = ^(3) (4) = 1, 5(3) (5) = 2, 5(3) (6) = 4, 5(3) (7) 



Remark 3.10 

The recurrences (3.14) do not have any polynomial solutions or solutions in hypergeometric terms for 
I > 0. This fact follows by applying the algorithm Hyper, described in j21j . 

Finally, we now return to our main theme: the log-convexity of secondary structure numbers. 

Theorem 3.11 

The sequence [S^^\n)) is log-convex. 
Proof 

We start from the short recursion (3.15). Dividing (3.15) by S'^^''(n— 1) and denoting S'^^''(n)/S'^^^(n— 1) 
by Qn, we obtain the following recursion for the numbers Qn- 



1 



qn 



n + 2 



n — 1 2n — 5 
2n + 1 + + 



n - 4 



Qn-l Qn-lQn-2 qn-lQn-2Qn-3 . 



n > 4, 



(3.18) 



with the initial conditions qi = q2 = 1, qs = 2. It is easy to check that < (?2 ^ • • • ^ 
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Assume, for the moment (and we shall prove it later on) that the sequence is convergent with 
limit q when n — > oo. By passing to limit in (3.18), we obtain the equation q'^ — 2q^ — — 2gf + 1 = 0, 
whose maximal positive solution is q = (p'^ = ip + 1, where ip = ^'^^ is the golden ratio. 

Define now the sequence = 2n+3 • clearly an increasing sequence and its limit is (p'^. We 

claim that (a^) is interlaced with our sequence (qn)- More precisely, we shall prove by induction that 

an<qn< On+l, (3.19) 

for 77, > 6. 

First we check directly the cases n = 6, 7, 8 and 9. Now take n > 9. Prom the induction hypothesis 
and (3.18), we have 

, „ ,n— 1 n— 1 n — 4 n — A 

(n + 2)g„ = 2n + 1 + + + 



7n-l qn-iqn-2 qn-iqn-2 qn-iqn-2qn-3 

= 2n + H \ \-{n-4)- 



qn—i qn—iqn—2 qn-iqn-2qn-3 

n—1 n — 1 , a„_3 — 1 

> 2n + l + + + (ra-4)- 



We would like the right hand side to be at least (n + 2)a„. But this is equivalent to 

(2n + l)anan-ian-2 + (n - l)a„_i(o„_2 + 1) + (n - 4)(a„_3 + 1) - (n + 2)a^an_ia„_2 > 

Inserting the formulae for a„'s, we get 

12(5 + V5)n^ - 2(241 + 121^/5)71^^ + 2(847 + 382^/5)n^ - 3(341 + 146V5)n + 126 + 54^5 

(2n - 3)(2n - l)(2n + l)(2n + 3)2 

The denominator is positive for all integers n > 2. Denote the numerator by L{n) and shift its 
argument for 6. The polynomial L{n + 6) has only positive coefficients, so it can not have a positive 
root. It then follows that L{n) can not have a root 7 > 6. So the left inequality is valid for all n > 6, 
and hence qn>. o-n- 

23 



> 0. 



To prove the other inequahty, note that the induction hypothesis imphes 



n — 1 n — 1 n — 1 , , an-2 — 1 
(n + 2)gn < 2n + 1 + + + + (n - 4)- 



The condition that the right hand side of this inequahty does not exceed (n + 2)a,i_|_i is equivalent to 

(2n + l)a„_ian_2«.n-3 + (n - l)a„„3(a„_2 + 1) + (n - 4)(a„_2 - 1) - (n + 2)a„+ia„_ia„„20n-3 < 0. 

Substituting the formulae for a„'s, we get 

^ (82 + 42^/5)^3 - (572 + 248\/5)n^ + (1103 + 474^/5)n - (529 + 247\/5) ^ ^ 
(2n - 3)(2n - l)(2n + l)(2n + 5) " 

If we put n + 5 instead of n in the numerator, we get a polynomial with all the coefficients positive, 
and from this we conclude that the numerator does not change the sign for n > 6. So, we have proved 
the inequality qn < a^+i, and thus completed the induction step. This proves the theorem. ■ 

Corollary 3.12 

The sequence qn = strictly increasing for all n > 5, bounded from above by 93^ and 

lim^^oo qn = 'p'^ = ■ 

Remark 3.13 

It is proved in [12j that the asymptotic behavior of S^^\n), n > 0, is given 

where Ki and ai are constants depending only on the rank /. Denote qn^ = S^^\n) / S^^\n — 1). Then 
it follows 

/ 1 \ 3/2 

qll'^aiil--] /ai. 

In other words q}} asymptotically behaves as an increasing sequence tending to a; as n — > cxd. It can 
be shown that ai G [2, 3]. The exact values of ai are known for I < 6. So, oq = 3, ai = (3 + \/5)/2, 
02 = 1 + \/2, etc., and a; \ 2 as Z ^ 00. 
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Theorem 3.14 

The sequences S^''^ (n) are log-convex, for / = 2, 3 and 4. 



Outline of the proof 

We present only the case I = 2. Dividing the short recursion for S^'^\n) by S^'^\n — 1) and denoting 
the quotient gfa) ^^["''-^.^ by qn \ we obtain the recursion for the sequence (qn^) 



'in 



n + 2 



n — 1 n — A 2n — 11 n — 7 

2n + l + + 



(2) ' (2) (2) (2) (2) (2) (2) (2) (2) (2) (2) (2) (2) 

'n-1 1n-l1n-2ln-3 1n-l1n-2ln~3^n-4 Qn-l1n-2'ln-3^n-4Qn-5 



(2) (2) (2) (2) (2) 

with the initial conditions Qi — q2 — Q3 = ^, Q4 = = We want to prove that the sequence 
{qn^} is increasing. 

Prom Remark 3.13 we conclude that the sequence (qn^) behaves asymptotically as (1 + -\/2)(l — ^)^^'^- 
Denote this quantity by 6„, i.e. 6„ = (1 + \/2)(l — ;^)^^^- Now take the first three terms of the series 
expansion of 6„ in powers of . Define 

3 1 3 1 \ 8n2 + 20n + ll 

"2 1 - — + -7 —ZTK = OL2- 



2 n + 2 8 (n + 2)2 y 8(n + 2)2 ' 

where a2 = 1 + \f2. The sequence (a„) tends increasingly toward a2- We shall show now that the 
sequences (q'i^'') and (a„) are interlaced, i.e. 

for sufficiently large n. 



(2) 

Suppose inductively that a„_i_i < < dn-i for i = 1, 2, 3, 4, 5. Then 



n + 2 



n—1 n — 4 2n — 11 n — 7 

2n + 1 + + 



an-2 an-20.n-3an-4 an-ian-20n-3an-4 an-ian-20n-3an-4an-5 



(2) 

If we prove that the right hand side of this inequality does not exceed a„, the right inequality, qn < an, 
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will follow. But this is equivalent to the condition 

Pioin) 



>0, 



Qi2{n) 

where Pio{n) and Quin) are certain polynomials in n of degree 10 and 12, respectively. (Using 
Mathematica, the polynomials Pio{n) and Qi2{'n) can easily be computed explicitly.) Their leading 
coefficients are 262144(H-\/2) and 8^(l + \/2)^, respectively, and we can conclude that this quotient is 
positive for all n big enough. Again, the biggest real roots of the polynomials Pio{n) and Qi2{'n) can 
be easily found using Mathematica. It turns out that their quotient becomes (and remains) positive 
for n > 39. 

(2) 

On the other hand, from the induction hypothesis and the recursion for g„ ' it follows that 



7P) > — 



n + 2 



„,n— 1 n — 4 2n— 11 n — 7 

2n + 1 + + 



fln-l a„_ia„_2a„_3 a„_2an-30n-4an-5 an-2an-30'n-40'n-50'n6 



That the right hand side of this inequality is > a„_i is equivalent to q^J^^) > 0, where P13 and Q15 
are certain polynomials in n of degree 13 and 15, respectively, with the positive leading coefficients. 
Their quotient is positive for n big enough (n > 6). 

The claim now follows by checking that {q^'^\n)) is increasing for n < 44. 

A similar proof works for Z = 3 and 1 = 4. We omit the details. ■ 



4 Calculus method 



Let us start again as in (3.1) with a linear homogeneous recursion for positive numbers a(n), and 
consider again the corresponding recurrence (3.2) for the quotients q{n) = a{n)/a{n — 1). Suppose 
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again we want to prove that a(n) is log-convex, i.e. that {q{n)) is an increasing sequence (at least 
from some place uq on). 

This time we do the following. Define a continuous function / : [1, oo) — R (or / : [no,oo) — R) 
starting with the appropriate linear function on the segment [1,2] (or [no, no + 1], for some no G N) 
determined by the initial and the next value of q{n), and then continue to the next segment by the 
same rule as (3.2). In other words, by replacing q f , n x, (3.2) becomes 

fix) = R{x) + -jp^y (4.1) 

defined so for x > 2 (or x > no + 1). In a sense, / is a dynamical system patching the discrete values 
/(n) = q{n) if R and S are "good" enough functions. For example, if R and S are rational functions 
(and in combinatorics this is mostly the case) without poles on the positive axis, then / is a piecewise 
rational function, i.e. rational on every open interval (n, n+ 1) for any integer n > 1 (or n > no). If we 
can prove that / is smooth on such open intervals (usually by proving that / is bounded from above 
and below by some well-behaved functions), we can consider the derivative f'{x), for any x G (n, n-|-l). 
The idea is to show inductively on n that f'{x) > (or f'{x) < if we want to prove log-concavity 
of a(n)). This will imply that / is an increasing function (or a decreasing function in the log-concave 
case) on any open interval (n, n -|- 1), and then, by continuity of / it will follow that / is increasing 
(or decreasing) on its whole domain, and hence q{n) = f{n), n G N, will increase (or decrease), too. 

Now, in general, if we want to prove that a sequence q{n) defined by (3.2) is increasing, we form the 
corresponding functional equation (4.1) with the appropriate start, i.e. we define f{x) on some starting 
segment [no, no -|- 1] to be an increasing function and then prove inductively on n that f'{x) > for 
X G (n, n -|- 1). It is always necessary to have some a priori bounds for /. So, suppose we know 
< m{x) < f{x) < M{x), for all x > uq. Let us find some sufficient conditions which ensure that 
f'{x) > 0. Of course, we assume that R and S are smooth on all open intervals (n, n -|- 1), n > no- 
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Fix an X G (n,n + 1) and write / = f{x), fi = f{x — 1), R = R{x), S = S{x). Then (4.1) can be 
written as 

h 

or equivalently 

ffi = Rfi + S. (4.2) 
Taking the derivative d/dx of both sides of (4.2), we get 

f'fi + ffi = R'fi + Rfi + S', 



implying 



From (4.2), we obtain 



f'h = R'h + S' + {R-f)fi. 



f'fi = R'fi + S'-^f[. (4.3) 



/: 



1 

Assume that 

< m{x) < f{x) < M{x), (4.4) 

for all X > no, and suppose inductively that / is increasing on [no, n]. This means that /{ > 0. Then, 
with our notation convention mi = m{x — 1) etc., we have: 

Theorem 4.1 

If R' > 0, R'mi + S" > and 5 < 0, then /' > 0. 
Proof 

Obvious from (4.3), (4.4) and the above discussion. ■ 
Theorem 4.2 

Suppose R' >0, S' >0, S >0 and mim2{R'mi + S') > S{R[M2 + S[). Then /' > 0. 
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Proof 

Divide (4.3) by /i and substitute in by the same rule /(. WE obtain 



, _ E!h + S' S , _ R'h + S' S 

and this imphes 



-^1/2 + S[ Si , 
/2 h 



fi _ R'h + S' S R[f2 + S[ SSi I 

By the inductive hypothesis, /2 = f'{x — 2) > 0, and since SSi > 0, it follows that the last term in 
(4.5) is non-negative. On the other hand, 

E^h + S^ _ S R[h + S[ ^ ^ ^ ^^^^^^^^ ^ ^ ^ 

Jl /l J2 

Since R' , S' , S > and mim2{R'mi + S') > S{R[M2 + S[), by our assumptions, it follows that (4.6) 
holds, and we are done. ■ 

Now let us show how to apply the above method to some combinatorially relevant numbers. 

Recall that the n-th big Schroder number rn is the number of lattice paths from (0, 0) to (n, n) with 
steps (1,0), (0, 1) and (1, 1) that never rise above the line y = x. Equivalently, r„ is the number of 
lattice paths from (0, 0) to (2n, 0) with steps (1, 1), (1, —1) and (2, 0) that never fall below the x-axis. 
For n = we put ro := 1. So, is the number of Motzkin paths whose every plateau, every valley, 
every terrace and every plain is of even length. Hence, via Proposition 3.6, big Schroder numbers 
count secondary structures whose unpaired bases form contiguous blocks of even length. The n-th 
little Schroder number is s„ = ^r„, sq = 1. 

Theorem 4.3 

The Schroder numbers are log-convex. The sequences rn/vn-i and Sn/sn-i increasingly tend to 3-|-2\/2 



as n — ^ 00. 
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Proof 

Similarly to the proof of Proposition 3.2, it can be shown (and it is well known) that (r„) satisfies the 
following (convolutional) recurrence: 

n 

j=0 

This implies that the generating function for (r„) is given by 



, , 1 — X — Vl — 6x + 
r(x) = . 

This, in turn, implies (as in Proposition 3.2) that r„ satisfy the following short recursion 

rn = [3(2n - l)r„_i - {n - 2)r„„2], (4.7) 

n + 1 

together with the initial conditions tq = 1, ri = 2. (There are also combinatorial proofs of (4.7), see 
or m].) Divide (4.7) by r^-i and denote r^/r^-i by q{n). We get 

, , 3(2n- 1) n-2 1 , , 

n + 1 n + 1 g(n — 1) 

with q{l) = 2. Then g(2) = 3 and q{3) = 11/3. Define the function / : [2,oo) ^ M by 

i(2x + 5) , xG[2,3] 

iiT[6^-3-7fE%] , x>3. 
It is easy to show by induction on n that / is bounded on [2, n] , for n > 2. More precisely, 3 < f{x) < 6, 

for all X >2. Also, / is continuous everywhere and differentiable on open intervals (n,n + 1), n > 2. 

In the notations of Theorem 4.1, we have 

R(x) = s{x) = -^^, m(x)=3, M(x) = 6. 

X + 1 X + 1 



fix) 



Then 
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Clearly, S{x) <0 for x > 2, and let us check that R'rrii + S' >0. But this is obvious, since 



9 
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(X + 1)2 (x + 1)2 (x + 1) 



> 0. 



Since f'(x) = 2/3 > for x G (2, 3), it follows from Theorem 4.1 that f'(x) > for all xe{n,n+l), 
n > 2. Hence, by continuity, / is increasing, so {q{n)) is increasing and therefore (r„) is log-convex. 
By passing to limit in (4.9) we get the last claim. ■ 

Now that we are more familiar with this method, let us show once again (but this time almost 
automatically) the log-convexity of Motzkin numbers. 

Recalling the short recursion (3.9), we define the function / : [2, oo) ^ M by 

2 , xe[2,3] 

^^[2x + 1 + 'f0^] , x>3. 



Here we have 



2x + l 3(x-l) 
R{x) = , b[x) — 



x + 2 



x + 2 ' 



so 



R'{x) 



> 0, S'{x) 



> 0. 



(x + 2)2 - ^ ^-^^ (x + 2)2 - 
Further, it is easy to check that 2 < /(x) < 7/2 for all x > 2. So, we may take m(x) = 2, M(x) = 7/2. 
Let us check the inequality mim2{R'mi+S') > S{R[M2+S[) from Theorem 4.2. But this is equivalent 
to 



2 • 



and this is equivalent to 



+ 



9 



(x + 2)2 (x + 2)2j - x + 2 



> 



3x — 3 



7 



+ 



9 



2(x + 2)2 (x + 2)2j ' 



1.5x2 + 61.5x + 117> 0. 

This last inequality is certainly true for all x > 1. Hence, by Theorem 4.2, /' > 0, and this implies 
the log-convexity of Motzkin numbers. ■ 
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A directed animal of size n is a subset 5 C (NU {0})^ of cardinality n with the fohowing property: 
if p £ S, then there is a lattice path from (0, 0) to p with steps (1,0), (0, 1), all of whose vertices lie in 
S. Let an be the number of directed animals of size n. It is well known that the generating function 
for (a„)„>i is given by 



A(x) = V a„x" = - ( A/ - 1 I =x + 2x'^ + 5x^ + 13x'^ + 35x^ + . . . 

These same numbers appear also as the row sums of Motzkin triangle (see ^U] and @]). 

Theorem 4.4 

The sequence (an)n>i of numbers of directed animals is log-convex. Hence the sequence an/an-i 
increasingly tends to 3. 

Proof 

Taking the derivative A'{x) of the generating function A{x) we get 

(1 + 2;)(1 - 'ix)A'{x) = 2A{x) + 1, 
and equating coefficients of x"^^ above yields the recurrence 

nan = 2na„_i + 3(n — 2)a„_2, > 3, 

with ai = 1, 02 = 2. 

The corresponding function for the successive quotients is given by 
and /(x) = ^ on [2,3]. It is easy to check that 2 < f{x) < |. 

Here R{x) = 2, S{x) = > and R'{x) = 0, S'{x) = ^ > 0, and m = 2, M = 7/2, so the claim 



follows from Theorem 4.2. 
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n 

'n\r 



We aplly now our approach to Franel's sequences. Recall that the n-th Franel number of order r is 
defined by 

sir) ^ 

k=0 

The numbers Sn^ satisfy a linear homogenous recurrence of order [^^J with polynomial coefficients. 
For r = 3, the recurrence reads as follows (HZ]): 



n 



= (7n2 -7n + 2)s2i + 8(n - l)^^^^, (4-9) 



with the initial conditions S^^ = 1, S^^ = 2. 



Theorem 4.5 

The sequence ( S^^M is log-convex. The same is true for the sequence ( S^^ 
\ ) n>0 V 



n>0 



Proof 

Again, taking the quotients q{n) = Sn / S\_t^ in (4.9), we get 

, , 7n2-7n + 2 8(n-l)2 1 ^ 
q{n) = 2 + 2 -( TT' "^2, 4.10 

with q[\) = 2. Then g(2) = 5. Define the function / : [1, oo) ^ M by 

3x-l , xG[1,2] 



/(^) 



7x^-7x+2 _L 8(x-l)^ 1 



^ 1 — j-^^ i\l 1 x>2. 



The function / is bounded. More precisely, 5 < f{x) < 9 for all x > 2. Indeed, / is clearly bounded 
on [1,2]. On [2,3] / is given by 

_ 21x2 -41x + 18 
^^^^ ~ 3x2 _ 4^ ' 

and one can check easily that 5 < /(x) < 9 on this interval. Now, by induction on n > 3, it is not 

hard to verify that if 5 < /(x) < 9 for x £ [2,n], then also 5 < /(x) < 9 for x € [n, n + 1]. Also, / is 

a continuous function. So, here we have, in the notations of Theorem 4.1: 

7x2 _ 7^ + 2 8(x-l)2 

R{x) = ^ , S{x) = K 

x^ x^ 
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and we may take m{x) = 5, M{x) = 9 for x > 2. Since R'{x) = > 0, S'{x) = > and 

S(x) > for X > 2, it remains only to check the last condition in Theorem 4.2, i.e. 
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7x-2 16(x-l) 

5 • 5 I 5 



> 



8{x - 1)' 



7x-9 16{x - 2) 



+ 



(X - 1)3 (x - 1) 



This is equivalent to 



643x2 - 1021x + 650 > 0, 



and this last inequality is certainly true for x > 1. 

So, as before, we conclude that / increases, hence (g(ra))^>Q is an increasing sequence and therefore 

( Sn^] is log-convex. 

V /n>0 

Note that {q{n)) converges, being increasing and bounded, and so, passing to limit in (4.10) we obtain 
that 



lim 



c(3) 



■'n-1 



A similar proof applies to the case of ( sii^ ] , starting with the 

V /n>0 



recurrence 



n 



S^^^ = 2[6n^ - + 5n - l<S\S^^l^ + (4n - 3)(4n - 4)(4n - b)S^^l^, 



(4) 



Sj^^ = 1, ^i^^ = 2. It turns out that 



We leave the details to the reader. 



lim 



c(4) 



s. 



16. 



n-l 



As a bit more simple case, let us show the log-behavior of sequences defined by two-term recurrences 
with constant coefficients. 
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Proposition 4.6 

Let {an)n>o ^ positive sequence defined by 

On = Cittn-l — C2a„_2, 

where Ci, C2 > are constants, and with some initial conditions oq > 0, ai > 0. Then the log- 
behavior of (on) is completely determined by the log-behavior of its first three terms. In other words, 
the sequence (a^) is log-convex if of < 0002, and log-concave if of > 0002- 

Proof 

Let q{n) = an/cLn-i- The recursion for q{n) is given by 

C2 



q{n) = Ci- 



g(n- 1)' 

and the first two values are ^(1) = ai/ag, q{2) = 02/01. Set A = 0002 — o-i and define the function 
/: [l,oo)^Rby 

m-- 

/ is a linear fractional mapping on any interval (n, -|- 1), n G N, and / is continuous. 



+ a , xG[1,2] 

aoai ^ ao ' L ' J 

C^ "FT — ^ 1 \ 1 X 2. 



Suppose first that A > 0. Then / is obviously increasing on [1,2] and bounded from below by 
/(I) = g(l) = ai/oQ. Namely, suppose that f{x) > /(I) on [l,n]. For x G [n,?7, + 1] we have 

and so f{x) is bounded from below on the whole interval [l,n -|- 1]. Hence / is bounded from below 
by /(I). This also shows that / is difFerentiable on any open interval (n, n -|- 1), n G N. Here we have 
R{x) = Ci, S{x) = —C2 (so R'{x) = S'{x) = for X N) and we may take m(x) = ai/uQ. It follows 
by Theorem 4.1 that /'(x) > for all x N. By continuity, / is an increasing function, and hence 
/(n) = q{n) is an increasing sequence. 
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If A < 0, then / decreases on [1, 2] and similarly as above / is, in this case, bounded from above by 
/(I). Since 

fix) = C2^^^^, 

we see that the sign of the derivative is propagated to the right and is determined by its sign on (1, 2). 
Hence, in this case, f'{x) < 0. ■ 

Corollary 4.7 

The sequence (i*2n+i)„>o of odd-indexed Fibonacci numbers is log-convex, while the sequence of even- 
indexed Fibonacci numbers (-F2n)„>i is log-concave. 

Proof 

Let an = i*2n+i- The numbers an satisfy the recursion a„ = 3a„_i — a„_2, and = 4 < 5 = 0902. The 
numbers bn = i*2n satisfy the same recursion bn = 36n-i — bn-2, but this time with initial conditions 
bi = 1, 62 = 3. However, 62 = 9 > 8 = bib^. So, the claim follows from Proposition 4.6. (Of course, 
this Corollary is an immediate consequence of the Cassini identity — = (—1)".) ■ 

As our next application of the method, we consider some classical orthogonal polynomials. Let ly > 
— 1/2 be a parameter. The Gegenbauer (or ultraspherical) polynomials Cn\t) are defined by 

m) 

nC^'\t) = 2t{u + n- l)ct\(t) - {2u + n - 2)C^^},{t), n>2, (4.11) 
together with the initial conditions C^\t) = 1, c'f\t) = 2ut. 

For V = 1/2, Cn^'^\t) = Pn{t) is the Legendre polynomial and for v = 1, Cn\t) = Un{t) is the 
Chebyshev polynomial of the second kind. 
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Theorem 4.8 

The sequence (Cn\t)j is log-concave loi t > 1, u > 1, and log-convex for < < 1, t > 

V / n>0 



max 



'n>0 

1 1 
/2u 



Proof 



^ ^*^(i-'-)-i ^+ ^*^(3'^-i)+i ,ifxG[l,2], 

(4.13) 



Let qt{n) = Dividing the recursion (4.11) by nC^\{t) we obtain the recursion 

, ^^v + n-12v + n-2 1 , 
n n ^t(n — 1) 

with initial condition gt(l) = 2vt. Then gt(2) = tiy + 1) — ^- Now define the function ft : [1, oo) — M 

by formula 

ft{x) 

The function /t(a;) is piecewise rational, i.e. it is rational on all segments [n, n -|- 1], n G N. It is easy 
to show, by induction on n, that ftix) is continuous at the points x = n, for all n G N, and that 
h{n) = qt{n). 

Let us first consider the case u > 1 and t > 1. We claim that ft{x) is smooth on all intervals {n,n + l), 
n G N. To prove this, it is enough to show that ft{x) is bounded from below by some positive quantity. 
Let us check that ft{x) > j for all a; > 1. 

Since u > 1, t > 1, the function ft{x) is non-increasing on the interval [1,2], so it is enough to show 
that ft{2) > J. But this is equivalent to t{i/ + l) > ^, or 2t^(i/+ 1) > 3. Since > 1, we have 4t^ > 3, 
and this is true for alH > 1. 

Suppose now inductively that ftix) > j for all x G [l,n], and take an x G [n,n + 1], n > 2. Prom 
(4.13) we have 

„ , , 2t{u + X-1) 2V + X-21 ^ 
ftix) > Y =t, 
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and t > i for all t > 1. So, the function ft{x) is nowhere zero on [l,oo), hence it has no poles on 

[1, oo) and hence is smooth on all intervals (ra, n + 1), ra G N. 

The function ft{x) is obviously decreasing on [1,2]. Suppose now that it is decreasing on [l,n]. Let 
X G {n,n + 1). Then 



The term in square brackets is positive, since ft{x) > j, so the whole first term in (4.14) is negative. 
The second term is negative by the induction hypothesis, hence //(.x) < 0. This completes the step of 
induction, and we can conclude that the function ft{x) is decreasing on [l,oo). Hence, the sequence 



induction on n that ft{x) > j on [l,oo), hence ft{x) is bounded from below by a positive quantity, 
hence ft{x) is smooth on all open intervals (n,n + 1), n G N. 

Suppose now that ft{x) is increasing on [l,n] and let a; € {n,n + 1), n > 2. Take a look at (4.14) 
again. The first term is positive since ft{x) > j, and the second term is positive by the induction 
hypothesis for all x >2. Hence fl{x) > and ft{x) is increasing on [1, n + 1]. The increasing behavior 



In Fig. 5 we can see the areas of log-convexity (denoted by A) and log-concavity (S) of sequences 




(4.14) 




is log-concave for all > 1, t > 1. 






Our analysis is not sufficient to determine the logarithmic behavior of sequences 




n>0 



whose parameters {u, t) fall into the areas denoted by I and H. 
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Figure 5: Areas of log-convexity (A) and log-concavity (S) for sequences of values of Gegenbauer 
polynomials 

Remark 4.9 

So far we have restricted our definitions of log-convexity and log-concavity to positive (or perhaps 
non-negative) sequences, since most combinatorial sequences have only non-negative terms. The 
definitions, however, make sense even if we include negative numbers. 

Corollary 4.10 

The sequence \Cn^(t)\ is log-concave for \t\ > 1, > 1, and log-convex for < i/ < 1, \t\ > 

V / n>0 

max < -k= , . ^ > . 
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Proof 

It follows from Theorem 4.7 and the relation 

cM(-t) = (-i)"cH(t), 

valid for all 7^ and all t. ■ 
Corollary 4.11 

The sequence (f^n(i))„>o of the values of Chebyshev polynomials of the second kind is log-concave for 
\t\ > 1. In particular, F2„, = Un{3/2) is log-concave. 

Proof 

Follows from Corollary 4.10, since Un{t) = ci^\t). ■ 
Corollary 4.12 

The sequence (-Pn(i))„>o of the values of Legendre polynomials is log-convex for \t\ > 1 and log-concave 
for Itl < 1. 



Proof 

Pn{t) = Cn^'^\t) and log-convexity follows from Corollary 4.10. The log-concavity follows from the 
well known fact that |-Pn(t)| < 1 for |t| < 1 and from the inequality (|^, 

valid for Itl < 1. ■ 



Our results on Legendre polynomials complement and generalize in some respects those of Turan and 
Szego (jSU)) who considered determinants of the form 



P„{x) P,^ + l{x) 

Pn-l{x) P„(x) 



Now, a combinatorial consequence is in order. The n-th central Delannoy number D{n) counts 
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the lattice paths in {x,y) coordinate plane from (0,0) to (n, n) with steps (1,0), (0,1) and (1,1). 
(Such paths are also known as "king's paths".) The generating function of (Z)(n))^>Q is given by 
D{x) = Qx+x ^ f \27\ ) . This implies easily that D(n)'s satisfy the recurrence (4.11) for u = 1/2 and 
t = 3 and -D(O) = 1, -D(l) = 3. Hence D(n) = -Pn(3), the value of the n-th Legendre polynomial Pn 
at t = 3. 

Corollary 4.13. 

The sequence D(n) of Delannoy numbers is log-convex and lim„^oo = 3 + 2y/2. 

Proof 

The log-convexity follows from Corollary 4.12. Hence qs^n) = D{n)/D(n — 1) is increasing. From 
(4.12) for 1^ = 1/2 and t = 3 it follows easily by induction that q3{n) < 6. Passing to limit, we obtain 
the second claim. ■ 

Remark 4.14 

No wonder that the above limit coincides with the limit for Schroder numbers from Theorem 4.3. This 
can be explained by showing that 

(For a proof, see e.g. PT].) 

The derivatives of ultraspherical polynomials behave log-concavely. More precisely, the following is 
valid. 

Theorem 4.15 

The sequence (Cn^ (t) ) is log-concave, for all > 0, t > 1. 

V / n>0 
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-^^{2ut^ + l)x + {4u+l)t + l , if [2,3], 

2tiu+x-l) 2U+X-1 1 If > q 

x-l x-1 ftix-1) , li X ^ O. 



Proof 

Deriving the recurrence (4.11) with respect to t and using the relation 
we get the two-term recursion for derivatives: 

which starts with C^''^' {t) = 0, c['''^' (t) = 2u, C^"^' (t) = Aiy{u + l)t and C^''^' {t) = 4u{u + l){u + 2)t'^ - 
2u{iy + 1). Hence, the successive quotients qt{n) = Cn^ {t)/C^\{t) satisfy the recurrence 

v + n-l 2i/ + n- l 1 

qtin) = 2t — , n > 3, 

n — 1 n — 1 qt{n — 1) 

starting with qt{2) = 2{y + l)t, q't(3) = (i^ + 2)t — ^. The appropriate function ft : [2, oo) — R is 
defined by 

Mx) = 

The function ft{x) is obviously decreasing on [2, 3] for all u > 0, t > 1. Also, ft{x) is continuous in all 
points a; = n, for n G N, and /t(n) = qt{n). As before, it is not hard to check that ft{x) > j, for all 
X >2. Hence, the derivative /[{x) exists for any x G (n, n + 1), ra > 2. 

This derivative is obviously non-positive on (2,3). Suppose inductively that //(x — 1) < 0. The first 
term above is negative since /t(a;) > j, and the second term is negative by the induction hypothesis. 
Hence, ft{x) is decreasing on (ra, n -|- 1), and then, by continuity, on [l,n -|- 1]. This completes the 
inductive step and hence gt(n) = /t(n) is also decreasing. ■ 

We conclude our review of orthogonal polynomials with Laguerre polynomials. The (ordinary) La- 
guerre polynomials satisfy the recursion 

nLn{t) = (2n - 1 - t)Ln-i{t) - (n - l)L„_2(i) 
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with the initial conditions Lo(i) = !> = 1 — i- 



Theorem 4.16 

The sequence (L„(t))^>Q is log-concave for all t <0. The ratio L„(t)/L„_i(t) is greater than one for 
alH < and it tends to 1 decreasingly as n ^ oo. 



Proof 



Dividing the above recursion by nL„_i(t) and denoting ^^^^ji^ by qt{n), we get a recursion for qt{n) 

2n-l-tn-l 1 -,r\ 

qt{n) = -T 4.15 

n n qt{n — 1) 

with qt{l) = 1 - t. By computing qt{2) = ^^^^i^, we see that qt{l) > qt{2), for all t < 0. Define 
ft : [1, oo) — R by formula 

2a:-l-i z-l 1 



ftix) 



7~7 TT , if ^ 2. 

X /t(x-l) ' — 

This function is, again, piecewise rational, continuous for x = n G N and /t(n) = gt(ra) for all n € N. 
By induction on n > 2, one can easily check that, for x G [1, n], we have 1 < ft{x) < 1 — t for alH < 0. 
Hence, 

l< ft{x) <l-t, x>l,t<0. (4.16) 
By computing the derivative, we see that //(x) < for a; G (1, 2), and for a; G (n, n+ 1), n > 2 we find 



x^ x'^ ft{x — 1) \ xj f1{x — l) 

Suppose f[{x — 1) < 0. Then for t < —1, all three terms on the right hand side are negative, hence 
fi{x) < 0. For — 1 < t < 0, the second and third terms are still negative, but the first term is positive. 
The claim will follow if we prove 

1+t _ 1 1 ^ ^ 
x'^ ft{x - 1) ~ 
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But this is equivalent to (1 + — 1) — 1 < 0, and this reduces to — 1) < ^pi . But ft{x) < l — t, 

and for —1 < t < we have 1 — t < j^. Hence, the sum of the first two terms is negative, and 
fl{x) < 0. By continuity, ft is decreasing, and so qt{n) is decreasing. Since by (4.16) this sequence is 
also bounded, it is convergent and passing to limit in (4.15) we get that this limit is 1 as claimed. ■ 

Some of our results on orthogonal polynomials, in particular those concerning Legendre and Laguerre 
polynomials, have been already known (see, e.g. JU), but here we derived them in a simple and 
unified manner, almost automatically. 

So far we have been considering only the two-term recurrences. The "calculus" method, however, 
works as well for higher order recurrences, as the following example shows. 

Let L{n) be the number of graphs on the vertex set [n], whose every component is a cycle, and define 
L(0) := 1. Then L{1) = 1, L(2) = 2, L(3) = 5, L(4) = 17 etc. The following recurrence holds (|23, 
Ex. 5.22): 

L{n) =nL{n-l) - (^^~^^ L{n-3), n > 3. (4.17) 

Theorem 4.17 

The sequence (L(n))^>Q is log-convex. 
Proof 

Let q{n) = L{n)/L{n — 1). From (4.17), dividing with L(n — 1), we obtain 

fn-l\ 1 

q{n)=n-[ — — — , n > 3, 

V 2 J q[n- l)q[n - 2) 

with the initial conditions q{l) = 1, q{2) = 2. The value g(3) is equal to 5/2 



(4.18) 

We claim that (g(n))„>i 
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is an increasing sequence. To this end, define the function / : [1, oo) — R by 

, ifxG [1,2], 

, ifxG [2,3], 



m = { 



x+2 
2 



(4.19) 



X — 



(x-l)(x-2) 



if X > 3. 



2 f(x-l)f{x-2) ' 

The function / is obviously continuous, rational on any [n,n + 1], n E N, without poles on [n, n + 1] 
and clearly increasing on [1,3]. Let us first prove an a priori bound 



f{x) > Vx + 1, x>2. 



(4.20) 



This is clearly true for x G [2,3]. Suppose that it is true for x G [2,n], and let a; G [n, n + 1]. Then 
from (4.20) we have 



{x-l){x-2) 1 (a;-2)V^^ 
j{x) > X _ , = X — 



2 y/x\/x — 1 



1^ 



and it is easy to check that 



(x - 2)VJ^ ^ , 

X - '-^ > a/xTT, X > 2. 

2Vic 

The function / is differentiable for all x > 1, x ^ N. Let x G (n, n + 1) for some n > 2 and denote 
= j{x - 1), /2 = /(x - 2). Suppose inductively that /{ > 0, > 0. Then 



2 /1/2 



The last term is positive, and from (4.21) it follows that 



1 _ ^ T~r ^ /1/2 ^ \fx\/x- 1 = \/x2 -X > X - ^, 

^ /1/2 ^ 



for X > 2. Hence, /'(x) > and theorem is proved. 



Consider now a general three-term recurrence of the form 

a(n) = R{n)a{n - 1) + S{n)a{n - 2) + T(n)a(ra - 3), n > 4. 
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(4.21) 



Dividing this recurrence by a{n — 1) and denoting q{n) = a{n)/a{n — 1), we obtain the recurrence for 

q{n): 

q{n) = R{n) + ^^+ n>3, (4.22) 

q{n — 1) q(n — l)q{n — 2) 

with some given initial conditions, say, g(l) = bi, q{2) = 62, 62 > ^'i > 0. 

Suppose we want to prove that {q{n))n>i is an increasing sequence. Again, form the function / : 
[1,00) — R mimicking the rule (4.22) and starting with the linear functions on [1,2] and [2,3], 
connecting the points (l,6i) and (2,62) and the points (2,62) and (8,63), respectively, where 63 = 
q{S) = i?(3) + + 1^ (supposing, of course, 63 > 62)- For x > 3, the function / is defined by 
replacing qhy f and n by x in (4.22). 

for X > 3 (or x > uq + 1, for some uq G N). For a fixed x we write simply /(x) = /, R{x) = R, 
S{x) = S, r(x) = T, f{x - 1) = /i, /(x - 2) = /2, etc. So, we write (4.23) simply as 

f = R+T + ^^ /1/2/ = Rhf2 + Sh + T. (4.24) 
/i /1/2 

We assume that i?, S and T are "good enough" functions, in the sense that / is differentiable on open 
intervals (n, n + 1) for integers n>2 (or n > no). This is usually a consequence of an a priori bound 
of the type 

< m(x) < /(x) < M(x), (4.25) 
for some well-behaved functions m(x) and M(x). 

Taking the derivative d/dx of both sides of (4.24) gives us 

^ (4.26) 



/1/2 



■ Sh + T T ■ 
^ 7 Ji - -7-J2 
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where 

F = R'fif2 + S'f2+T'. (4.27) 

Suppose that R,S,T > and R' , S' , T' > 0. By substituting the analogous expressions for /{ and 
/2 in (4.26) and supposing inductively that /g, and /4 > we obtain (denoting Fi = F{x — 1), 
Ri = R{x — 1), etc.) that: 

/1/2/' = F- ^Pj^Fi - y^F2 + Bf^ + C/^ + Dfi 
/1/2/3 J2J3J4 

where B, C and D are non-negative. To prove that /' > 0, assuming inductively that /g, f!j^ and 
/4 > 0, it is enough to prove the following: 

/1/2/3 72/3/4 ^ ~ 

or equivalently 

fMsUF > {Sf2 + T)UFi + r/iF2. (4.28) 
Taking into account (4.25), it suffices to prove a stronger inequality 

m^F > {MS + T)MFi + MTF2 (4.29) 

for all X > no, for some no G N, and to check that f'{x) > for x G (fc. A; + 1), A; < no- 

Let us now consider a concrete example. A Baxter permutation is defined in 21.. The numbers 
B{n) of Baxter permutations in S„ satisfy the recurrence 

(n + l)(n + 2)(n + 3)(3n-2)B(n) = 2(n + l)(9n^ + Sn^ - 4n + 4)5(n - 1) 

+ (3n - l)(n - 2)(15n2 - 5n - 14)5(n - 2) 

+8(3n + l)(n - 2)2(n - 3)S(n - 3), n > 4, (4.30) 

together with the initial conditions S(0) = 1, B{1) = 1, B{2) = 2, B(3) = 6. 
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Theorem 4.18 

The numbers B{n) of Baxter permutations are log-convex. The limit lim„_>oo 5(?7.)/i?(n — 1) exists 
and is equal to 8. 

Proof 

Let q{n) = B{n)/B{n — 1). Prom (4.30) we form the recurrence for q'(n)'s, and then according to the 
initial values q{l) = 1, q{2) = 2, q(3) = 3, we form the function / : [l,oo) — >■ M defined by f{x) = x 
on [1, 3] and for x > 3 by the rule 

(x + l)(a; + 2)(a; + 3)(3x-2)/(x) = 2(a: + l)(9a;^ + 3x^ - 4x + 4) 

(3a; -!)(.;:- 2)(15;r- - 5,r - 14) 

_^8(3a; - l)(x - 2f{x - 3) 



f{x-l)f{x-2) 

This, written in the form of (4.23), yields to conclude that R, S, T are positive rational functions with 
no poles on [l,oo). We also see that R{x) y 6, >S'(a;) y 15, T{x) 8 as x — > oo. 

It is apparent from (4.31) that f{x) is a piecewise rational function, i.e. rational on intervals [n, n + 1], 

n G N. For example, for x G [3, 4], 

_ IS-T^ + 51.T^ - 122.T''^ - 87x2 _^ 200x + 12 
~ (x- l)(x + l)(x + 2)(x + 3)(3x-2) 

Clearly, /(x) is continuous everywhere. It can be checked, by induction on n, that the function / is 
bounded. More precisely, 

7 < f{x) < 9, X > 47. (4.32) 

We want to prove that / is an increasing function. Prom the above a priori bound, it follows that / 
is differentiable on all open intervals (n, n + 1), n G N. The non-negativity of /'(x) can be checked for 
X < 49, using e.g. Mathematica, and for x > 49, it follows from the stronger inequality, 

7^F > (9 • 15 + 8) • 9Fi + 9 • 8F2, 
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obtained by substituting appropriate values for m, M in (4.29). As this inequality is true for all a; > 9, 
the first claim follows. The second claim follows passing to the limit in (4.30). ■ 



Many other three-or-higher-term recurrences can be investigated by this method. Let us only mention 
that the four-term recurrence (3.15) for the numbers S^^\n) of secondary structures of rank 1 can 
also be shown to be log-convex by this "calculus" method. The details (rather tedious) are given in 
jllj . jl3j . Another example is the number Sn oi n x n symmetric matrices with entries 0, 1, 2, whose 
sums of all rows and all columns are equal to 2. These numbers satisfy (HZj) the following recurrence 

Sn = (2n - l)5„_i - (n - l)(n - 2)5„_2 - (n - l)(n - 2)5„_3 + l(n - l)(n - 2)(n - 3)S„_4, 

starting with Sq = 1, 5i = 1, S2 = 2>, = 11. The sequence ('S'n)„>o is also log-convex. 

Finally, as it should be clear by now, this method applies to any P-recursive sequence (a(n)), satisfying 
a recurrence of the form 

Q{n)a{n) = Pd{n)a{n - 1) + Pd-i{n)a{n - 2) + . . . + Po(n)a(n -d-l), n>d+l, (4.33) 

where d > is an integer and Pq) • • • j Pd^ Q real polynomials, Q > 0. The corresponding function for 
the successive quotients q{n) = ^^^^i'^ is given by the functional equation 

Q{x)f{x)=Pd{x) + ^f^ + ... + - -, x>d+l. (4.34) 

f{x-l) f{x-l)...f{x-d) 

As we have seen, the most important thing in this approach is to express the derivative f'{x) in terms 

of previous derivatives. So, fix a point x £ {n,n + 1), n > d, and write as before for short / = /(x), 

fj = f{x — j), j = 1,2, . . ., Pi = Pi{x), Qi = Qi{x), i = 0,1, . . . ,d. Then (4.34) can be written as 

d p 

Qf = ~r~f f ' 

Jlh ■ ■ ■ Jd-i 

or by denoting the product of all values by 11, i.e. 11 = /1/2 ■ ■ ■ fd, and by 11^ = /1/2 ■ ■ ■ fj the partial 
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products (so 11(^ = 11 and Hq = 1), as 



UQf = Y,ifd~i+i---fd)Pi- 



(4.35) 



Taking the derivative d/dx of both sides in (4.35), after some manipulations, we obtain the following 
formula (in terms of Wronskians) : 

1 



^ " 1 1 A PiU' 



i=0 



Q P. 
Q' P[ 



^^^d-i 



In particular, for d = 1 this reduces to 



1 


Q 


Po 


1 

^ g2 


Q 


Pi 




Q' 


p' 


Q' 


Pi 



Qto ■ 



f/ 
72JI- 



(4.36) 



(4.37) 



With a priori bounds < m{x) < f{x) < M{x) and by substituting /( in (4.37), one gets almost 
instant proofs of log-behavior. For example, if Pq ^ and if we want to prove the log-convexity, 
hence assuming f[ > 0, then if the first Wronskian Wq in (4.37) is positive, we only have to check 
jjWq + Wi > 0, and check that / increases at the beginning. 

Of course, not every (combinatorially relevant) sequence satisfying a recurrence of this type can be 
expected to have a reasonable log-behavior; it is enough to recall here the sequences ek{n) from Section 
2, whose log-behavior is rather chaotic for A; > 3. 

Let us only mention here that the log-convexity of secondary structure numbers of general rank I 
can also be proved by calculus method, using the explicit formulae from Proposition 3.9 and formula 
(4.36). The details will appear elsewhere. 



As a finall remark, note that our approach applies also to linear nonhomogeneous recurrences for 
positive numbers. So, for example, let (a(n)) be given by the linear recurrence of the first order 



a(n) = R{n)a{n — 1) -|- <S'(ra) 
50 



(4.38) 



Consider the quotients g(n) = ^"^["''-^^ and note that 



a{n) = q{n)q{n-l)...q{2)a{l), n>2. 



(4.39) 



Then, dividing (4.38) by a{n — 1) we obtain a (long) recurrence for (;(n)'s: 



q{n) = R{n) + 



S{n) 



q{n-l)q{n-2)...q{2)a{iy 



(4.40) 



To get a short recurrence for q{nys, substitute for a(n) and a(n — 1) the corresponding products (4.39) 
in (4.38) (n > 3): 



qin)q{n - 1) . . . g(2)a(l) = R{n)q{n - 1) . . . g(2)a(l) + S{n) = 



^^^^ g(n)g(n-l)...g(2)Q(l) ^ ^^^^ (?(n)g(n - 1) . . . g(2)a(l) 



g(n) 



g(n)g(n-l)...g(2)a(l)' 



Prom there we get 



and then 



q{n)qin-l)...q{2)a{l) S{n) 



q{n - 1) . . . q{2)a{l) S{n-1) 



1 - 



_ Rin) 
lin) J ' 

R{n - 1) 
q{n - 1) J 



Substituting (4.41) in (4.40) yields a short recursion for g(n)'s: 



(4.41) 



q{n) = R{n) + 



5(n) R{n - l)S{n) 1 



S(n - 1) S{n - 1) q{n - 1) 



Similarly, for a second order linear recurrence 



a(n) = R{n)a{n - 1) + 5(71)0(71 - 2) + T{n), 



we obtain 



g(7z) = R{n) + , + ^ ^ 



i?(7l - 1) 5(71 - 1) 



q{n-l) T{n- 1)1 q{n - 1) q{n - l)q{n - 2) _ 



Then we can proceed as before. 
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5 Calculus method in two variables 



We shall outline our method for non-negative sequences a(n, k) in two integer variables n,k>0 (or 
n > no, k > ko). Suppose (as often in combinatorics) that the numbers a{n,k) satisfy a two-term 
recurrence of the form 

a{n, k) = R{n, k)a{n — 1, A; — 1) -|- S{n, k)a{n — 1, k), (5.1) 

with some known functions R and S, together with some initial values, usually of the type a(0, 0) = a, 
a(l,0) = b, a(l,l) = c. Suppose we want to prove that the sequence a{n,k) is log-concave in k, i.e. 
that a(n, k)'^ > a{n, k—l)a{n, k + 1), for all n, k. Here is what we do. Write down (5.1) with k replaced 
hy k-1: 

a{n, k-l) = R{n, k - l)a{n -l,k-2) + S{n, k - l)a{n -l,k- 1). (5.2) 

Denote 

a{n,k) 
a{n, k — 1) 

and divide (5.1) by (5.2) (always assuming we do not divide by zero). 

R{n, k)a{n -l,k-l)+ S{n, k)a{n - 1, k) R{n. k) + S{n, k)q{n - 1, k) 



q{n, k) 



R{n, k - l)a(n - 1, A: - 2) + S{n, k - l)a{n - l,k - I) ^^"^7^ + Sin, k - I) 

q(n—l,k—l) ' V ' / 

R{n,k) + S{n,k)q{n-l,k) 

= q{n — l,k — Ij- 



i?(n, k-l) + S{n, k - l)q{n -l,k-l) 
Equivalently, 

g(n, k)[R{n, k - 1) + S{n, k - l)q{n -l,k-l)]= q{n -l,k- l)[R{n, k) + S{n, k)q{n - 1, k)\. (5.4) 

The log-concavity of a(ra, A;)'s is equivalent to q{n,k) > q{n,k + 1), for any fixed n and all k. The 
idea is again to pass to a "continuation" of (5.4) by letting n^x, k^y, q^f and obtaining the 
functional equation 

fix, y) [R{x, y-l) + S{x, y-l)f{x-l,y- 1)] =f{x-l,y- l)[R{x, y) + S{x, y)f{x -l,y)]. (5.5) 
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We assume that R and S are "good enough" functions, in the sense that / is continuous everywhere 
and smooth on open cells (n, n + 1) x (m,m + 1), for all m,n. What we want to prove is that / is 
decreasing in y for any fixed x. Fix a point (x, y) in an open cell Q = {n,n + 1) x (m, m + 1), and 
prove inductively that 



5/ 
dy 



x,y) < 0. 



(5.6) 



For the fixed pair {x,y) write for short fij = f{x — i,y — j), for i,j = 0, 1,2, . . ., and similarly for R 

and S. So, /oo = f{x,y), Roi = Rix,y — 1) etc. In this notation, (5.5) can be written as 



f[Roi + Soifu] = fii[R + Sf: 
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(5.7) 



Now take the partial derivative d/dy of both sides in (5.7). We have 



^[i?oi + 5oi/ii] + /[^ 



f dSoi „ dfn 
JU—^ r OQi 



dy 



df 



11 



df 

[Rqi + 5'oi/ii] = /ii 

dy 



dR^^ dS_ 
dy dy_ 



dy _ 

dRio 



dy 



dy 

dSoi 



[R + Sho] + fii 



OR ^ _^^^dS ^ ^dfw 
dy dy dy 



11- 



dy 



f + ^fnS + ^[R + Sho - Soif] 



dy 



dy 



Substituting here / from (5.7), we get 



5/ 
dy 



[Roi + Soifii] 



f: 



11 



Roi + Soifu 



\Ro^ + ^oi/ii)(f + /lof ) - (R + ^/io)(^ + /ii^ 



H 



dy 



Roi + Soifu dy 



(5.8) 



Assume that R and S are positive. Hence if / starts with some positive values, then / can be considered 
positive, too. Suppose inductively that 



dfio dfu 
dy ' dy 



< 0. 



Then the last two terms in (5.8) are negative and to prove (5.6), it is enough to prove that the first 
term is negative, too. In other words, to conclude (inductively) that (5.6) holds, it is enough to prove 
that the "free" term is non-positive, i.e. that 



(-Roi + Soifu) ( — + ho-Q^ 



(R + Sf: 



\ dy dy 



< 0. 



(5.9) 
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So, if we can check that / begins decreasingly in y and assuming inductively that / is decreasing in 
y, then by (5.9) we can conclude that / is decreasing in y at the point and then, by continuity, 

that / is decreasing in y everywhere. 



Note that the following inequalities imply (5.9) (simply by comparing similar terms): 

^ dR dRoi ^ OR dSoi ^ dS dR^i ^ OS dSoi 
ay oy ay oy ay ay ay oy 

In terms of Wronskians, writing G' for these inequalities can be written in the form: 



Ri 



01 



R'n 



01 



R 
R' 



< 0, 



01 



01 



R 
R' 



< 0, 



Roi S 
^01 ^' 



<0, 



5*01 S 

or 0/ 
'-'01 



< 0. 



(5.10) 



Instead of formalizing everything (which can be done with a little care), let us take an example. 
Example 5.1 

The Eulerian number E{n,k), is the number of permutations vr from E„ with exactly k ascents, 
i.e. with exactly k places where irj < ttj+i- We know (see, e.g. ^Hj) that these numbers satisfy the 



recurrence 



E{n, k) = {n- k)E{n -l,k-l) + {k + l)E{n - 1, k), 
with the initial conditions E{0,k) = Jofcj E(n,0) = 1, n,/c > 0. Then (5.4) becomes 

q{n, k)[n — k + 1 + kq{n — 1, A; — 1)] = q{n — 1, k — l)[n — k + {k + l)q{n — 1, k)] 

with q{n, k) = E(n, k)/E(n, k — 1), for 1 < k < n. The initial conditions are q{2, 1) = 1, g(3, 1) = 4, 
g(3, 2) = 1/4 and (we extra define) q{2, k) = 0, for k > 2. We want to prove that q{n, k) > q{n, k + 1), 
for any fixed n and all k. In the sense of the above discussion and notations, here we have R{n, k) = 
n — k, S{n, k) = k + 1. Passing to the natural "continuation", (5.7) becomes 



f[x-y + l + yfn] = fii[x - y + (y + l)/io]. 



(5.11) 
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In fact, we define the function / : [2, oo) x [1, oo) — M first on two shaded strips in Fig. 6 below and 
then continue by the rule (5.11). In Fig. 6 we indicated the values /(n, k) = q{n, k) ioi n > 2, k > 1 
in the lattice nodes. On the vertical walls of the (square) cells Q2, Q3, Q4, ■ ■ ■, as well as on their 
lower horizontal walls, we define / to be appropriate linear functions. The upper walls oi Q3, Q4, . . ., 
are determined by (5.11). The right walls of Q3, Q4, . . . are also determined by (5.11). We fill in / on 
cells Qn, Q'n by appropriate homotopies connecting (possibly nonlinear) functions on the walls. 
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Figure 6: The boundary values of function / 



So, let /(x, 2) = 2"(a; - n) + 2" - - 1, for x G [n,n + 1], n > 2, /(2, y) = 2 - y for y € [1, 2] and 
/(2, y) = for y > 2. Further, let /(3, y) = for y G [1, 2], and, in general, let f{n, y), y G [1, 2], 

be the linear function in y between the points (n, q{n, 1)) and (n, q{n, 2)). The cell Q2 is surrounded 
by linear functions o(x) = 3x - 5, x G [2,3], 6(y) = ^i^, y G [1,2], c(x) = x G [2,3] and 
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d{y) = 2 — y, y e [1, 2]. Extend / to Q2 by the homotopy 

fix, y) = (2 - y)a{x) + (y - l)c(x) = ^{23x + 18y - llxy - 38), (x, y) G Q2. 

The cell Q'3 is surrounded by functions a'^ix) = ^,x e[2,3], h'^{y) = 43^:^^, V e [2, 3], ^(x) = 0, 
X G [2,3] and (i3(y) = 0, y G [2,3]. Extend / to by the homotopy 

/(x,y) = (3 - x)4(y) + (x - 2)6:,(y) = ^^'^^^'fj^ , {x,y) G g^,. 

Extend / to Q'^, Q'^, ... to be zero. Next, extend / to Q^, Q4, ... by the homotopies connecting the 
lower walls of Qn, given by a„(x) = 2"(x — n) + 2" — n — 1, for x G [n, n + 1], n > 2, and the upper 
walls, given by rational functions Cn{x) determined inductively on n by (5.11), thus obtaining / by 

Qn 

f{x, y) = (2 - y)a„(x) + (y - l)c„(x), (x, y) G Qn- 

For example, since a2(x) = 3x — 5, a3(x) = 7x — 17, C2(x) = then 

a2(x - l)[x - y + (y + l)c2(x - 1)] _ 3x — 8 



X — y + 1 + ya2(x — 1) 4 



and hence / 



is given by 



. , , . 53x + 60y-25xy-128 
/(x,y) = (2 - y)a3(x) + (y - l)c3(x) = . 



In this way, / is well defined on the shaded strips on Fig. 6 and extended to [2, cxd) x [1,oo) by the 
rule (5.11). 

It is easy to check that / is continuous and nonnegative and that it is a rational function with no 
poles on any open cell, and hence smooth on any open cell. It is also easy to check inductively on n 
that / is decreasing in y on and Q'^ for n > 2, i.e. ^{x,y) < for (x,y) G int{Qn) U int{Q'^). 

Now that we have elaborated carefully the "beginning" of /, the rest is more-or-less automatic. The 
inequality (5.9) reduces to (x - y + 1 + y/ii)(-l + /lo) < (x - y + y/10 + /io)(-l + /ii)- This is 
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equivalent to 



1 + xfio + fio < -fio + xfii + /lo/ll 



or, after some rearrangement 



/ii/io - 2/10 + 1 + x(/ii - /lo) > 0. 



The second term is non-negative by the induction hypothesis, and the rest is non-negative since 



Hence, -^{x, y) < for (x, y) G int{Q), for ah Q. So, / is decreasing in y on every open cell, and hence 
by continuity, / is decreasing in y everywhere. In particular, g(n, k) = f{n, k) > f{n, k+1) = q{n, k+1) 
and we are done. 

In the same manner we can prove that ^{x — y + 1) > /, hence ^ > 0, for any {x,y) £ int{Q), for 
any Q, and this implies q{n + l,k) > q{n, k), for any fixed k and all n. Hence, 



Apart from settling the "beginning" of /, we can (almost automatically now) prove the well-known 
log-concave behavior in the second variable of the binomial coefficients (and in general find the log- 
concave behavior when R and S in (5.1) are constants), g-binomial coefficients, Stirling numbers of 
the first and second kind, Schlafli numbers, cover many particular results (e.g PUI) and so on. 

Of course, the method can be extended in a few ways; for example, to three-or-more term recurrences, 
recurrences for three or more variables a{n, k, I) etc. but we shall not consider it here. 



/ii/io - 2/10 + 1 > f?o - 2/10 + 1 = (/lo - 1)' > 0. 




E{n + l,k) E{n + l,k-l) 
E{n,k) E{n,k-1) 



> 0. 



In a word, a general idea of this method is as follows. Combinatorics gives a recurrence. Pass to 
the quotients of the neighboring members (in the variable under consideration), pass to the natural 
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"continuation" /, find some bounds (upper, lower or both, depending on the nature of the problem) of 
/, make sure that / is differentiable on open cells, examine the rate of growth of / at the "beginning" 
(i.e. check the sign of the derivative there), and finally, prove inductively from the associated functional 
equation corresponding to the recurrence, that the sign of the derivative remains the same. 



58 



References 



[1] M. Aigner, Motzkin numbers, Europ. J. Combinatorics, 19(1998)663-675. 

[2] M. Abramowitz, I. Stegun, Handbook of Mathematical Functions with Formulas, Graphs and 
Mathematical Tables, Appl. Math. Ser. 55, U.S. Govt. Printing Office, Washington, 1964. 

[3] E. Babson and E. Steingrimsson, Generalized permutation patterns and a classification of the 
Mahonian statistics, Seminaire Lotharingien de Combinatoire, B44b(2000)18pp. 

[4] E. Barcucci, R. Pinzani, R. Sprugnoh, The Motzkin family, PU.M.A. Ser. A Vol. 2 (1991) 
249-279. 

[5] E.A. Bender and E.R. Canfield, Log-concavity and related properties of the cycle index polyno- 
mials, J. Comh. Theory A, 74(1996)56-70. 

[6] M. Bona, R. Ehrenborg A combinatorial proof of the log-concavity of the numbers of permutations 
with k runs, J. Comh. Theory A, 90(2000) 293-303. 

[7] F. Brenti, Unimodal, log-concave and Polya frequency sequences in combinatorics, American 
Math. Society, Providence, RI, 1989. 

[8] D. Callan, Notes on Motzkin and Schroder numbers, preprint, 2000. 

[9] A. Claesson, Generalised pattern avoidance, preprint, 2000. 

[10] R. Donaghey, L.W. Shapiro, Motzkin numbers, J. Comb. Theory A, 23(1977) 291-301. 

[11] T. Doslic, Problems of Matching Enumeration and some AppHcations to Biochemical Graphs, 
Ph. D. Thesis, University of Zagreb, Zagreb, 2001. 

[12] T. Doslic, D. Svrtan and D. Veljan, Secondary structures, submitted, 2001. 

59 



[13] T. Doslic and D. Veljan, Calculus proofs of some combinatorial inequalities, submitted, 2001. 

[14] M. Eweida, On Turan's determinant for Legendre and Laguerre polynomials. Rev. Math. Hisp. 
Amer. 15(1955) 79-87. 

[15] D. Foata, D. Zeilberger, A classical proof of a recurrence for a very classical sequence, J. Comb. 
Theory A 80(1997) 380-384. 

[16] R.L. Graham, D.E. Knuth and O. Patashnik, Concrete Mathematics, Addison- Wesley, Reading, 
1988. 

[17] C. Haslinger and P.F.Stadler, RNA structures with pseudo-knots: Graph-theoretical, combina- 
torial and statistical properties. Bull. Math. Biology, 61(1999)437-467. 

[18] D.E. Knuth, The Art of Computer Programming vol. 1, Addison- Wesley, New York, 1973. 

[19] J. Kruskal and D. Sankoff, Time Warps, String Edits and Macromolecules (2nd edition), Addison- 
Wesley, Reading, 1999. 

[20] D.C. Kurtz, A note on concavity properties of triangular arrays of numbers, J. Comb. Theory A 
13(1972) 135-159. 

[21] M. Petkovsek, Hypergeometric solutions of linear recurrences with polynomial coefficients, J. 
Symbolic Computation 11(1998) 1-22. 

[22] B. Sagan, Inductive and injective proofs of log-concavity results Discr. Math., 68(1988)281-292. 

[23] P. Schuster, I.L. Hofacker and P.F. Stadler, Combinatorics of RNA secondary structures, Discr. 
Appl. Math., 88(1998)207-237. 

[24] L.W. Shapiro and R.A. Sulanke, Bijections for Schroder numbers. Mathematics Magazine, 
73(2000)369-376. 

60 



[25] R. Stanley, Log-concave and unimodal sequences in algebra, combinatorics and geometry. Ann. 
N.Y. Acad. Sci., 576(1989)500-535. 

[26] R.P. Stanley, Enumerative Combinatorics vol. 1, Cambridge Univ. Press, Cambridge, 1987. 

[27] R.P. Stanley, Enumerative Combinatorics vol. 2, Cambridge Univ. Press, Cambridge, 1999. 

[28] R. Stanley, Positivity problems and conjectures in algebraic combinatorics, in Mathematics: 
Frontiers and Perspectives. (Eds. V. Arnold et al.), IMU-AMS, 2000, 295-319. 

[29] P.R. Stein and M.S. Waterman, On some new sequences generalizing the Catalan and Motzkin 
numbers. Discrete Mathematics, 26(1979)261-272. 

[30] O. Szasz, Inequalities concerning ultraspherical polynomials and Bessel functions, Proc. Amer. 
Math. Soc (1950) 256-267. 

[31] G. Szego, Orthogonal Polynomials, AMS, New York, 1959. 

[32] G. Szego, On an inequality of P. Turan concerning Legendre polynomials. Bull. Amer. Math. Soc 
54(1948) 401-405. 

[33] M.S. Waterman, Secondary structures of single stranded nucleic acids. In G.C. Rota, editor. 
Studies on Foundations and Combinatorics. Advances in Mathematics Supplementary Studies, 
Vol. I, Academic Press, New York, 1978, 167-212. 

[34] M.S. Waterman, Applications of Combinatorics to Molecular Biology, In R. Graham, M. Grotschel 
and L. Lovasz, editors. Handbook of Combinatorics, Elsevier, Amsterdam, 1995. 



61 



